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We argue that the degree of understanding the causes and mechanisms of the long-term variations 
(11-year and 22-year) in the galactic cosmic rays (GCR) characteristics is still insufficient and 
to improve it we need new approaches and methods. For the time being there is a long-lasting 
controversy on how these long-term variations, observed for more than 50 years in the inner 
heliosphere, are formed. It is widely believed that the 11-year variation is due entirely to the 
toroidal branch of solar activity (the area and number of sunspots, the strength of the heliospheric 
magnetic field etc) because of the diffusion, convection and adiabatic energy loss, while the much 
smaller 22-year variation is caused by the particle drifts connected with the poloidal branch of 
solar activity (the high-latitude solar magnetic fields). At the same time, both past and more 
recent numerical simulations indicate that the contribution of particle drifts could be significant 
for both 22- and 11-year variations in the GCR intensities. However, even those who agree on 
the significant influence of drifts appear to have different perceptions on the mechanisms of this 
influence. 

In this paper, we present an analysis of the possible causes of the first point of view (small role 
of drifts in the 11-year GCR variation) and the reasons why one can expect the significant contri¬ 
bution of the processes connected with the poloidal branch of solar activity in both types of the 
long-term variations of the GCR characteristics. Then we briefly discuss some numerical methods 
suggested in the past and recently and the approaches and perspectives for the sought-for methods 
are considered. 
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1. Introduction 


The long-term (or more speeifically the 11-year and 22-year) variations of the GCR intensity and 
anisotropy, connected with the toroidal (or sunspot) and the poloidal (high-latitude solar magnetic 
fields) branches of solar activity, have been studied for more than fifty years. Many of the proper¬ 
ties of GCRs are known from observations mostly in the inner heliosphere. The transport partial 
differential equation (TPE) serving as a theoretical basis for understanding and modeling these 
variations was formulated about fifty years ago [p 


21 , 


, which was followed by extensive 
work on heliospheric models and in-depth analysis of the TPE coefficients, describing the main 
processes (diffusion, convection, adiabatic energy loss or gain, and particle drift) involved in the 
modulation of GCR intensity. The systematic efforts of theoretical and numerical modeling have 
resulted in a reasonable overall description of the observational data [^, 16], although there are 


some questionable features in the simulations [19]. 


However, the degree of understanding the causes and mechanisms of the long-term variations 
in the GCR characteristics is still insufficient. As a cause of these variations we mean the change 
of some heliospheric parameters (such as solar wind velocity, strength of the magnetic field, ifs 
polarity efc.), connected wifh one of fwo branches of solar acfivify or bofh of fhem, which gives 
rise fo fhe GCR variafion under discussion. If should be sfressed fhaf fo call fhe above change of 
fhe heliospheric parameters fhe cause of fhe GCR variafion if is nof sufficienf fo demonsfrafe fhe 
regression befween fhem, buf if is necessary fo poinf ouf and undersfand how fhe mechanism of 
such influence acfs (see, e.g., ]|^). In sfudying fhe GCR variations fheorefically, we imply fhaf 
fheir mechanisms are confained in fhe boundary-value problem for fhe GCR intensify or, probably 


more practical, in fhe stochastic differenfial equations (SDEs) pSQ and in fhe coefficienfs of fhese 
equations (fhe diffusion fensor, fhe solar wind and drifl velocities and so on), used fo model fhe 
intensity . However, when one solves fhis fheorefical basic problems - fhe boundary-value problem 
or SDEs - to describe fhe observed long-term variations of fhe GCR intensity, one does nof yef 
obfain a full answer on fhe relafive imporfance of fhe differenf mechanisms and how if changes in 
differenl times, and in differenf regions of energy and space. To obfain fhese answers one needs 
some additional means and approaches for freafing fhe fheorefical basis for fhe GCR modulafion. 

Eor fhe time being fhere is a long-lasfing confroversy on fhe main causes and mechanisms of 
fhe long-ferm variafions of fhe GCR characferisfics. There is a widespread view (especially among 
fhe observers), fhaf fhe 11-year variation, J\]_{t), is almosf enfirely due fo fhe sunspof cycle (if 
is often even called fhe sunspof cycle, Jss{t), in GCR infensify), while fhe much smaller 22-year 
variafion, J 22 {t), is connecfed wifh particle drifls. As fo fhe GCR anisofropy, ifs 22-year varying 
componenf, W-wave, was fhe firsf observational facf fhaf fhe GCRs are sensitive fo fhe polarify 
of fhe high lafifude solar magnefic fields And fhe firsf driff model was devoted fo 

descripfion of fhis W-wave. However, fhe second componenf of fhe anisofropy, 11-year varying 
V-wave, was fhoughf fo be due enfirely fo fhe sunspof cycle. 

On fhe ofher hand, fhe active role of drifls in forming fhe 11-year variation of fhe GCR infensify 
followed from several model calculations (mosfly using firsf generation drifl models) [ ]I^ pT| , [T5| , 
28, The confribufion of differenl mechanisms fo fhe solar modulation of fhe GCR intensity 
was also scrutinized in several recenl papers (see [rs] , ^ ^^). Besides, recenlly a new 

numerical fool was suggesfed fo disenfangle fhe effecls of all fhe individual mechanisms from fhe 
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full calculated intensity [17 


In this paper we first consider the reasons for the different points of view on the causes and 
mechanisms of the long-term variations of the GCR intensity and anisotropy, then discuss some 
suggestions of the approaches and methods for improving the quantitative understanding of these 
variations. The possible perspectives for the sought-for methods are also considered. 


2. On different points of view 

There are several causes of the widespread belief that the 11-year variations in the GCR inten¬ 
sity and anisotropy are almost entirely due to the sunspot cycle The main cause is that (1) there 
is an overall anticorrelation between the GCR characteristics and different solar and heliospheric 
indices changing in phase with the solar activity cycle (the sunspot area, the strength of the regular 
heliospheric magnetic field (HMF) and ifs inhomogeneifies, efc.). Besides if is imporfanf fhaf (2) 
in fhe firsf 20 years (1955-1975) if was believed fhaf fhere was only fhe 11-year sunspof cycle in 
fhe heliospheric parameters and fhere was nof any dominafing HMF polarify A (fhe sign fhe reg¬ 
ular HMF radial componenf in fhe N-hemisphere). After thaf 722(0 considered as small 
variation wifh respecf fo 7ii(f) as (3) fhe Earfh is nof in a favorable posifion fo observe fhe drifl 
effecls [^. Furfhermore, if is offen presumed (wrongly) fhaf (4) fhe deciding cause of a powerful 
11 -year cycle of GCR flux musf be a powerful sunspof cycle alone while fhe change of fhe polar¬ 
ify of fhe HMF, fogefher wifh fhe accompanying change of fhe heliospheric currenf sheef (HCS) 
befween low and high fill can only resulfs in small changes in fhe infensify of GCRs. Finally, if is 
also offen presumed fhaf (5) fhe 11- and 22-year cycles in fhe GCR infensify should be caused by 
fhe causes of fhe same periods. We should admif fhaf fhe significanf role of driff is still quesfioned 
from observations even during fhe minima of solar cycles [|3|] (see fhe discussion in ||^]). 

The alfemafive view, fhaf bofh 11-year- and 22-year variafions in fhe GCR characferisfics are 
caused by fhe combined acfions of fhe sunspof cycle in fhe heliospheric characferisfics and change 
of fhe HMF polarify fogefher wifh fhe HCS filf, fhe laffer being imporfanf for bofh GCR variafions, 
is based on fhe assurance fhaf (1) fhe change of fhe dominafing HMF polarify A wifh 22-year cycle 
is a global process in fhe whole heliosphere and (2) fhere are mechanisms of fhe 11-year variafion 
in fhe GCR intensify forming as fhe 2-nd harmonic connected wifh fhe cause, changing wifh fhe 22- 
year period. Among fhese mechanisms we can mention: i) modulation of fhe 22-year mechanism 
wifh an 11-year period. For example, fhe drifl velocity depends on fhe HMF polarify (changing 
wifh 22-year period), and on fhe HMF slrength and fhe form of fhe HCS (bofh changing wifh 11- 
year period); ii) simultaneous change of two factors in the 22-year mechanism. For example, the 
direction of the drift velocity and the main channel of particle‘s arrival are changing synchronously; 
iii) two 22-year mechanisms acting in tandem. For example, if the GCR intensity is first modulated 
by some hypothetical mechanism external to the part of heliosphere inside the termination shock 
and then by drifts in this part, acting in phase with each other (both as sin(t (>22 • f)), then some part 
of the modulated intensity changes as sin^(ft )22 •?) = (! — cos(ft)n •t))/2, where (O 22 = lullly, 
coii = 27r/lly. 
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3. Some numerical tools and approaches 


To discuss the numerical methods for studying the role of different mechanisms in the long¬ 
term variations of the GCR characteristics we should first formulate the boundary-value prob¬ 
lem they deal with, first introduced in [25, 21, 26]. Usually instead of GCR intensity J{f,T,t) 
the boundary-value problem and TPE are formulated for the phase-space distribution function 
f{r,p,l) = J{f,Tf)/p^, where p is the momentum of particles. For the stationary case the TPE 
balances the divergences of the diffusion, drift and convection fluxes in space and thaf due fo adi¬ 
abatic cooling in the momentum space. For the case of axial symmetry in a spherical system of 
coordinates the boundary-value problem, i.e., the transport equation with boundary and "initial" 
conditions (which may also be written somewhat differently), appears as follows: 


■ ^ = _V.(^.V/) + v^w.v/-^—+y^'-.v/ = 0, 


diffusion 


convection-i-adiabatic loss 


drift 


(3.1) 


d r 


- 0 , f\r=r^—fnm{p), 


= 0 


)>=0,7r 


f\p=— fnm {Pmax ), 


(3.2) 

(3.3) 


where V^'^, V'^'' and are the solar wind and magnetic drift velocities, and symmetric part of 
the diffusion tensor, respectively. Besides, fnm{p) is the distribution function corresponding to 
nonmodulated GCR intensity and r^ax, Pmax are the radii of the inner and outer boundaries of 
the modulation region and the momentum above which there is no modulation. 

Besides the GCR intensity we are also interested in the long-term variations in the GCR 
anisotropy ^ = 38/{Anp'^vf), connected with the differential streaming S. The latter can be easily 
calculated after solving the boundary-value problem Eqs. (^-^), S = AKp^ (CcafV^'^ — K ■'V , 


where Ccg = ~l/35(ln/)/5(lnp), K and v are the Compton-Getting factor, the full diffusion ten¬ 
sor, and the particle’s velocity, respectively. In this paper we shall concentrate on the influence of 
the particle drift on the GCR intensity, as the role of drifts in forming the 11-year V-component in 
the anisotropy depends on their role in the 11-year variation of the GCR intensity. 

The importance of the particle drifts for the GCR intensity was proposed in 0 soon after 
the HMF model appeared of two unipolar magnetic “hemispheres” divided by the thin global HCS 
& The models [|T^, |I|, |T^, Q were the first ones getting both the 11-year- and 22-year long-term 
GCR variations from changing only the HCS tilt. These models had very specific predictions, for 
instance, marked differences in the GCR intensity time-profiles in consecutive solar cycles (plateau 
vs peaked shape in A > 0 and A < 0 cycles for positively charged particles, respectively), as well 
as differences in the radial and latitudinal gradients of GCRs. These features were confirmed by 
observations ii- Their magnitude, however, tend to be smaller than prediction obtained by 
purely drift-dominated models Some suppression of the drifts and modification of the TPE 
coefficients (the compound approach, 0) helps in describing many of the observed features in the 
GCR behavior quantitatively. 
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Already the first full drift calculations [10, ^ showed that drift tended to reduce modulation 
by significant factor for both polarity cases and now it is a common knowledge for the modeling 
community (see [1^, 24]). Recently in [1^], first, the parameters of the models were chosen in 
such a way that the calculations basically described the latitude, radial and energy dependencies 
of the observations for the HMF polarities, respectively, in 1987 and 1997 solar minima. Then the 
calculations were repeated with the same parameters but without the drift, that is omitting the drift 
term in Eq. (^T|). The calculated intensity almost everywhere decreased by a factor of 3-5 when 
compared with the intensities for both A > 0 and A < 0 periods. In [ jT^ it was shown that for the 
GCR intensity near the Earth the same was true for all other phases of solar cycle except the years 
near solar maxima. So the conclusion was made that in the models used, similar to [|^, ^ p^ ], 
the drift contribution was large for both polarities and it was emphasizes that particle drifts result 
in the significant part of the 11-year variation of the GCR intensity in the whole heliosphere. 

In [|^ devoted to modeling the unusual 23/24 minimum between solar cycles 23 and 24 (2009) 
the extent to which diffusion and particle drifts contributed to the total observed modulation from 
2006 to 2009 was also investigated. Eor this purpose the 3D stationary TPE ( |3.1| ) was solved by the 
finite-difference method. Besides the run with changes of all modulating parameters, the authors 
made the multiple runs each time changing only one of the factors, (e.g., the HCS tilt), while every 
other modulation parameter was kept unchanged. We shall not discuss the conclusions made by 
the authors, as here we are interested only in their approach to find the contribution of individual 
processes. In [ |^ also dealing with the unusual 23/24 solar minimum approximately the same 
method (the multiple runs each time changing only one of the factors) was used also to find the 
relative importance of different modulating factors. The main difference from the method used in 
[30] is that in [36] the time-backward Markov stochastic process method, proposed in [34, 
was used to solve the 3D SDEs. The questions of how the drift and other processes influence the 
propagation times and energy losses in the heliosphere were addressed in [|3^ also using SDE 
approach. The same method was used in [^ also changing one factor at a time (switching on and 
off the drift term and modifying the diffusion tensor) in multiple runs. 

In [ |r^ the method was suggested to decompose the calculated GCR intensity into the partial 
intensities connected with the main physical processes (the diffusion, convection, adiabatic cooling 
and magnetic drift) of the solar modulation, J = + Jp’"'’ + ■ In more details 

the method is described in [^], where it is also demonstrated how it can be used for studying 
the mechanisms of the GCR intensity modulation and the heliospheric structure during the solar 
minima with opposite HME polarities (A > 0, A < 0) and without drift effects (A = 0). 


4. Discussion and conclusions 

It looks like there is a hint of consensus between the modelers on the significance of particle drifts 
for both 11-year and 22-year variations in the GCR intensity. So when one discusses the observed 
drift effects, not only the difference between the GCR characteristics for A > 0 and A < 0 periods 
should be taken into account, but also the probable significant contribution of drifts in these charac¬ 
teristics for both polarities. As it is impossible to switch off particle drift in nature, estimating this 
contribution from the observations is a difficult task. Probably, drifts should also strongly influence 
the 11-year V-wave in the GCR anisotropy. 
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However, what we need are the numerical means for the quantitative study how the differ¬ 
ent modulation mechanisms form the long-term variations in the GCR intensity. Note that in all 
models [^, studying the contribution of different mechanisms by the method of multiple 


runs changing only one of the factors at a time, the overall distribution of the intensity gradients 
for each such run is not the same as it is in the run with all the modulating factors changing in 
time simultaneously. This difference in the intensity gradients can be very significant and so the 
diffusion, drift and convection of the particles, that is all the processes depending on the gradients, 
carry out differently in each of the runs. On the other hand, the decomposition of the calculated 
GCR intensity into the partial intensities connected with the main physical processes in |jn, 18, 20] 
is carried out in the run with all the modulating factors changing in time simultaneously, not dis¬ 
turbing the overall distribution of the gradients and hence the modulating processes. Probably, it 
is one of the causes of different results on the relative contribution of the drift and diffusion to the 
GCR intensity near the Earth during the unusual solar minimum 23/24 which were obtained in [ |^ 
and [|r^. According to ||3^ the relative contribution of the drift and diffusion is approximately 
the same, while in [ [T^ the drift partial intensity is always much smaller than the diffusion partial 
intensity the latter being almost compensated by the large negative partial intensity connected with 
the adiabatic energy losses. It should be said that the results of |]J] strongly depend on the models 
used (too simple and outdated in [|T^, |T^ |T7|, and there are unresolved questions in the method 
of partial intensities However, what we are looking for is the numerical instrument, which 

can be applied to any model. 

As up to now there is no accord on the ways how the long-term variations in the GCR char¬ 
acteristics are formed, even those who agree on the significant influence of drifts have different 
perceptions on the mechanisms of this influence. The first picture is that the GCR modulation is 
directly associated with the energy loss, and the latter depends on the path and dwelling time of 
GCRs in the heliosphere, these quantities being determined by the interplay of the processes in¬ 
volved. The drifts reduce modulation for both A > 0 and A < 0 periods, since drifts assist GCR 
particles to penetrate into the inner heliosphere faster in both polarity cases (either along the pole or 
at the HCS) ||T|, ||, Another picture arising from the consideration of the energy dependence 
of the gradients and the partial GCR intensities in the whole heliosphere [^j, is that the changes 
in the intensity distribution start at the highest energy in the outer heliosphere and accumulate with 
decreasing energy and radial distance. The main result is the change of the intensity gradients 
which influences the diffusion and other fluxes in such a way that the GCR intensity increases for 
both HMF polarities. 

The consensus in both pictures is that drift becomes more important and more apparent at 
larger helioradii and higher energies, and the global modulation is largely decided in the remote 
heliosphere (beyond tens of AU). Besides the main role of the energy loss in the first picture agrees 
with the conclusion in [l^j that the largest negative partial intensity is 7“^“^ (at least for cases 
considered there). It is also probable that the main difference in two above pictures is that in the 
first case we think of individual particles, while in the second picture we discuss the formation of 
the particle distribution: the two differs, but at the end could come with the same result. 

We also hope that both the finite-difference and stochastic methods can help. Moreover, prob¬ 
ably the method to improve the quantitative understanding of the long-term variations in the GCR 
characteristics consists in the synthesis of the two approaches. 
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